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By paying attention to the hole-doped two-dimensional 
systems of antiferromagnetically (strongly) correlated elec- 
trons, we discuss the cause of hole-rich phase formation in 
association with phase separation. We show that the phase di- 
agram obtained from the Maxwell's construction in the plane 
of temperature vs. hole density is consistent with one derived 
from the evaluation of hole-rich and electron-rich phases in 
real space. We observe that the formation of a hole-rich phase 
is attributed to the aggregation of hole pairs induced by spin 
singlet pairs present in the pseudogap phase and that a di- 
rect involvement of correlations between hole pairs are not 
essential for phase separation. 

PACS numbers: 64.75.+g, 71.10.Hf, 73.22.Gk, 74.80.-g 

Currently there has been an activated interest in the 
study of stripe modulations owing to strong experimen- 
tal evidences observed [|T|— |l3|| involving hole-doped high 
T c cuprates. Phase separation, a phenomenon of bro- 
ken symmetry at low energies can be manifested as a 
stripe phase or a phase-separated inhomogeneity. At a 
critical temperature, the phase separation appears ow- 
ing to the divergence of compressibility or the zero in- 
verse compressibility. The onset of phase separation can 
be obtained from the evaluation of the zero inverse com- 
pressibility using the well-known Maxwell's construction. 
Thus far, various numerical approaches ]l4|-|32| to the 
Hubbard and t — J Hamiltonians have successfully re- 
vealed the existence of phase separation in the hole-doped 
two dimensional systems of strongly correlated electrons. 
They include the Hatree-Fock (HF), the Lanczos exact 
diagonalization (LED), the density matrix renormaliza- 
tion group (DMRG), a fixed node quantum Monte Carlo 
(FNQMC) and Green's function Monte Carlo (GFMC) 
methods. Most recently a slave-boson functional integral 
(SBFI) method of Gimm and Salk has been proposed 
to reveal a good agreement with earlier exact numerical 
studies [|l5],|5) of phase separation boundary. However it 
remains to be seen whether the phase separation obtained 
from the Maxwell's construction can, indeed, show phase- 
separated distributions of charge and spin in real space. 
Currently there exists a great controversy ]3l| , |3^ , f2l"| , p5[ 
over the issue of whether the phase separation accom- 
panied by the hole-rich phase arises as a consequence 
of correlations between pairs of holes or not. Here we 
address this issue, based on our slave-boson functional 



integral approach |3(],|33|], and take the advantage of its 
usefulness to study the role of charge and spin degrees of 
freedom on phase separation. 

The t — J Hamiltonian to deal with the hole doped 
systems of antiferromagnetically correlated electrons is 
written, 

=-*E (4c* + H.c.) 
~ o E (4r c ij ~ 4 c 1t) ( c ^ c *t - c iT c a) • ( 2 ) 

In order to explicitly see physics involved with the spin 
and charge degrees of freedom, we rewrite the above 
equation in the f (1) slave-boson representation with the 
use of Ci a — fiab\ as a composite of spinon and holon 
hi operators subject to the single occupancy constraint 

H t -j = "* E flfjb]bi - | E kb 3 b)b\ 

The above U(l) slave-boson expression can also be de- 
duced from our earlier SU(2) slave-boson representation 
HH of the t — J Hamiltonian. The second term shows 
spinon pairing (spin singlet) interactions between adja- 
cent sites. Both terms explicitly manifest coupling be- 
tween the charge (holon) and spin (spinon) degrees of 
freedom. 

We attempt a multi-faceted study of phase separation 
by employing various levels of approaches; the Hatree- 
Fock (HF) |33| of Hubbard Hamiltonian for a qualitative, 
comparative study of phase separations, the Monte Carlo 
diagonalization (MCD) |34] for an accurate study of cor- 
related electron systems and the slave-boson functional 
integral (SBFI) method for an investigation of the role 
of charge (hole) and spin degrees of freedom on phase 
separation. Extending our earlier HF approach p3| of 
Hubbard Hamiltonian to a study of temperature depen- 
dence of charge and spin distributions in real space [ p5| , 
we first examine the domain of phase separation in the 
plane of temperature vs. hole density and compare the 
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phase separation domain derived from the Maxwell's con- 
struction with the one obtained from the direct (real) 
space calculations of electron- rich and hole-rich phases. 
Based on HF calculations with U = At, in Fig. |] we dis- 
play the predicted phase separation diagram in the plane 
of the temperature T vs. hole density x with the two 
different approaches; one from the use of the Maxwell's 
construction (denoted by shaded area in Fig. |l|) and the 
other from the direct evaluations of hole (charge) and 
electron (spin) distributions (denoted as solid diamonds) 
in the real space of 14 x 14 lattice. Temperature de- 
pendent HF total energies for Maxwell's construction are 
obtained from the solution of the following self-consistent 
equations (Eq. (j^)) involving the quasi-particle energy, 
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with e k = \J [-2t(cos k x + cos k y )} 2 + (^f) and 

|(c^Ci| — c^c^)|, by determining the uniform staggered 
magnetization m and the chemical potential fj, 
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Encouragingly the two results showed a large overlapping 
domain of phase separation in the T — x plane. The 
present mean field results are at best qualitative. 

In order to show the doping dependence of phase sep- 
aration in real space, we now examine the electron and 
hole distributions interior and exterior to the phase sep- 
aration boundary. At a chosen hole density of x = 0.25 
near but below a critical value x c ~ 0.3 at T = O.Oli, 
a stripe phase no longer appears and, instead, an inho- 
mogeneous phase separation is revealed, as is shown in 
Fig. 0. Although not shown here, such inhomogeneous 
structures are prevalent near but below the critical hole 
density x c , but they eventually disappear to yield uni- 
form phase structures beyond x c . It should be noted 
that despite the fact that no correlations exist between 
the Hatree-Fock quasiparticles, the phase separation is 
seen to occur as well known. They are not accurate for 
strongly correlated (large U) electron systems. However, 
it is to be noted that with currently available, numerically 
exact methods, the present level of real space calculations 
of temperature dependent phase separation for a large 
square lattice are not easily feasible. It is of great inter- 
est to see whether there exists a possibility of 'metastable' 
phase-separation in a certain region of temperature and 
doping where density fluctuations are excessively large 
as is shown by the results from the Maxwell's construc- 
tion. This case happened near T — O.Oli and x — 0.25, 
showing a metastable region of phase separation. In this 



small domain, the density fluctuations were predicted to 
be large. A detailed study for verification is necessary in 
the future. 

For an accurate account of strongly correlated elec- 
trons for the study of phase separation as a function 
of antiferromagnetic coupling strength J, we now take 
the MCD method Q of the t - J Hamiltonian, Eq. 
(1) with J = 1M for a 4 X 4 square lattice. Fig. || 
shows the computed ground state energy (in units of t) 
shifted by a linear factor e#n e as a function of the elec- 
tron density n e (n e = 1 — x). The solid triangle repre- 
sents the GFMC calculation by Hellberg and Manousakis 
p5| c and the open circles, our present results. From 
the Maxwell's construction the critical electron density 
(hole density) is predicted to be n c = 0.684 (x c — 0.316) 
compared to the value of n c — 0.730 (x c = 0.270) ob- 
tained from the GFMC. The overall variation of curva- 
ture as a function of electron (hole) density is grossly sim- 
ilar between the two approaches, although the predicted 
ground state energies are not the same. In Fig. |[ we 
display our MCD predicted phase separation boundary 
in the plane of electron (hole) density vs. antiferromag- 
netic interaction strength (Heisenberg coupling constant) 
in units of t and make comparison with other methods. 
The MCD agreed very well with the GFMC of Hellberg 
and Manousakis [§5| c , the LED of Emery et al. |l5| and 
the SBFI (U(l) slave-boson functional integral) result of 
Gimm and Salk (3^] . A salient feature is that all of these 
methods yielded a similar phase separation boundary, 
by showing a smoothly decreasing (increasing) trend of 
the critical electron (hole) density with an intercept near 
J/t ~ 3.5 as the antiferromagnetic interaction strength 
J ft increases. 

However, the accurate MCD calculations can not read- 
ily resolve the current controversy over the issue of 
whether the effect of correlations between the hole pairs is 
the primary cause of forming the hole-rich phase. For this 
cause we now explore the U(l) slave-boson functional in- 
tegral approach of the t — J Hamiltonian (Eq. (|J)). This 
method is advantageous to examine the role of the charge 
and spin degrees of freedom or the role of the holon and 
spinon degrees of freedom on phase separation. Earlier 
we reported that this method |30| also showed a satisfac- 
tory phase separation boundary for all J values (J/t < 1 
and J/t > 1) in general agreement with other numer- 
ical studies @|||. The spin (spinon) degrees of free- 
dom shown in the second term in Eq. (|^) allows spinon 
pairing (spin singlet formation) interactions between ad- 
jacent sites. Indeed it has been shown from the U(l) 
slave-boson theories that the spinon pairing, that is, the 
spin singlet pair order appears below the pseudogap (spin 
gap) temperature T* 
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33J . This indicates that the mo- 
tion of paired holes rather than the independent motion 
of separated holes is energetically preferred in the pres- 
ence of surrounding electron spin singlet pairs (spinon 
pairs) below the pseudogap (spin gap) temperature T*. 
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This is because the hole pairs (holon pairs) can readily 
migrate to the occupied sites of spin pairs (spinon pairs) 
involving no spin-bond breaking. 

In summary, we showed that the phase separation dia- 
gram obtained from Maxwell's construction in the plane 
of temperature vs. hole density is consistent with one de- 
rived from the real space (direct lattice) calculations of 
hole-rich and electron-rich phases; inhomogcneous phase 
separation appears near x c at finite temperatures; the 
formation of the hole-rich phase for phase separation is 
attributed to the aggregation of hole pairs induced by 
spin singlet pairs which exist in the spin-gap phase. 
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FIGURE CAPTIONS 

FIG. 1 The temperature (in the unit of hopping energy 
t) dependence of phase separation as a function of 
hole doping density x based on the Hatree-Fock 
approximation with U = At for a 14 x 14 square 
lattice. The shaded region represents phase sepa- 
ration obtained from the Maxwell's construction; 
the solid diamond, the phase separation obtained 
from the numerical evaluation of charge and spin 
distributions in real space and the open diamond, 
the uniform phase. 

FIG. 2 Spin (electron) and charge (hole) distribution at 
x = 50/196(~ 0.25) T = O.Oli. The large circles 
represent holes (charges) and smaller ones, elec- 
trons (spins). The size of circles represents the de- 
gree of hole-richness or electron-richness. 

FIG. 3 Maxwell's construction (dotted line) from predicted 
ground state energies denoted by open circles; the 
energies (in units of t) are shifted by a linear factor 
—enn,,, (eh(x) — enn e ) /t as a function of electron 
density for J = l.Ot with a 4 x 4 lattice. For com- 
parison the Green's function Monte Carlo result 
p5|] d (denoted by HM) with various lattice sizes 
up to 11 x 11 is displayed. 

FIG. 4 Phase separation for the hole-doped systems of an- 
tiferromagnetically correlated electrons in the plane 
of Heisenberg coupling strength, J/t and the elec- 
tron density, n e = 1 — x. The triangles denoted 
by HM represent the Green's function Monte Carlo 
prediction of Hellberg and Manousakis J2?| c and 
the stars denoted by EKL, the exact diagonaliza- 
tion result of Emery et al. p5[ . The solid cir- 
cles represent the results from the U(l) slave-boson 
functional integral approach of Gimm and Salk |3C| l 
and diamonds, our present MCD results. 



4 



0.3 



0.2 



0.1 



<> o o o ♦ 

<> ♦ 

♦ ♦ iiilllllV 



<>♦♦♦♦♦ 



0.1 0.2 0.3 0.4 

X 

FIG. 1. 



© © 


© 


© 


CD 


CD CD © 


© 


© 


© 


CD 


CD 


CD 


© CD 


CD 


© 


© 


CD © © 


CD 


CD 


© 


© 


CD 


© 


CD <J> 




CD 


© 


© © CD 


<t> 


* 


CD 


© 


© 


© 


CD <t> 


* 


CD 


© 


© © CD 


* 




CD 


© 





© 


© CD 


CD 


© 


© 


CD © © 


CD 


CD 


© 


© 


CD 


© 


© © 


© 


© 


CD 


CD CD © 


© 


© 


© 


CD 


CD 


CD 


CD © 




CD 


CD 


<f> CD CD 


© 


© 


CD 


CD 


<t> 


CD 


© © 


© 


© 


CD 


CD CD © 


© 


© 


© 


CD 


CD 


CD 


© CD 


CD 


© 


© 


CD © © 


CD 


CD 





© 


CD 


© 


CD <t> 


<t> 


CD 


© 


© © CD 


<t> 


<t> 


CD 


© 


O 


© 


CD <J> 


<t> 


CD 


© 


© © CD 


<t> 


* 


CD 


© 


© 


© 


© CD 


CD 


© 


© 


CD © © 


CD 


CD 


© 


© 


CD 


© 


© © 


© 


© 


CD 


CD CD © 


© 


© 


© 


CD 


CD 


CD 


CD © 


© 


CD 


CD 


<t> <D CD 


© 


© 


CD 


CD 


4> 


CD 



EKL Phase Separated 




FIG. 2. 



FIG. 4. 



0.25 




-0.25 



-0.75 



0.2 0.4 0.6 0.8 1 
electron density n e 



FIG. 3. 



5 



